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ABSTRACT 

The co-planarity of solar-system planets led Kant to suggest that they formed from an accretion disk, and 
the discovery of hundreds of such disks around young stars as well as hundreds of co-planar planetary systems 
by the Kepler satellite demonstrate that this formation mechanism is extremely widespread. Many moons in 
the solar system, such as the Galilean moons of Jupiter, also formed out of the accretion disks that coalesced 
into the giant planets. We report here the discovery of an intermediate system OGLE-2013-BLG-0723LB/Bb 
composed of a Venus-mass planet orbiting a brown dwarf, which may be viewed either as a scaled down version 
of a planet plus star or as a scaled up version of a moon plus planet orbiting a star. The latter analogy can be 
further extended since they orbit in the potential of a larger, stellar body. For ice-rock companions formed in 
the outer parts of accretion disks, like Uranus and Callisto, the scaled masses and separations of the three types 
of systems are similar, leading us to suggest that formation processes of companions within accretion disks 
around stars, brown dwarfs, and planets are similar. 

Subject headings: planetary systems - brown dwarfs - gravitational lensing: micro 


1. INTRODUCTION 

Brown dwarfs are intermediate in mass between stars 
and planets, making them important laboratories for test¬ 
ing theories about both classes of objects, including theo¬ 
ries about atmospheres, binarity, and formation mechanism. 
Brown dwarfs are also potentially a laboratory to test the¬ 
ories of planet formation, but detecting planets that orbit 
brown dwarfs is challenging. Four super-Jupiter planets or¬ 
biting brown dwarfs hav e been detected to dat e by three tech¬ 
nique s: direct-imaging (IChauvin et al.l l2004t iTodorov et al.l 
l2010l) . rad ial-velocity (iJoergens & Mulled l2007h . and mi- 
crolensing (iHan et al.l2013h . In the first three cases, the super- 
Jupiter planet was within a decade in mass of its host and well 
separated from it, suggesting that these systems more resem¬ 
ble scaled down binary stars than they do the extreme mass- 
ratio systems that are likely to emerge from accretion disks. 
In the fourth case, the separation was substantially closer and 
the mass ratio somewhat more pronounced so that the planet 
could have formed through the core-accretion scenario that is 
thought to be the origin of the solar system’s gas giants. 

With present technology, substantially lower-mass compan¬ 
ions to brown dwarfs, composed primarily of ice and rock, 
can only be discovered using the gravitational microlensing 
technique. Microlensing does not rely on light from either the 
host or planet, but rather infers the existence and properties of 
these bodies from their deflection and magnification of light 
from a more distant star that is fortuitously aligned with the 
system. As such, the low-luminosity of the brown dwarf and 
almost complete absence of light from its low-mass compan¬ 
ion do not interfere in any way with the discovery process. 

In this paper, we report the microlensing discovery of an in¬ 
termediate system OGLE-2013-BLG-0723LB/Bb composed 


of a Venus-mass planet orbiting a brown dwarf. The system 
may be viewed either as a scaled down version of a planet 
plus star or as a scaled up version of a moon plus planet orbit¬ 
ing a star, suggesting that formation processes of companions 
within accretion disks around stars, brown dwarfs, and planets 
may be similar. 

2. OBSERVATION 

Microlensing event OGLE-2013-BLG-0723 was discov¬ 
ered by th e Optical Gravitat ional Lensing Experiment 
(OGLE-IV) (lUdalski et akllMBl) on May 12, 2013 (HID ~ 
2456424.5) using its 1.3m Warsaw telescope at the Las Cam- 
panas Observatory in Chile. At (a,5)j2ooo=(17^34™40't51, 
-27°26'53".2) [{l,b) = (-0°.02,2°.83)], the event lies in one 
of nine fields toward the Galactic bulge that OGLE observes 
> lhr“*, which is sufficient cadence to discover teiTestrial 
planets. In fact, the discovery was triggered by a short ^ 
1 day, ^ 40% brightening of the lensing event that was still 
in progress at the time of the announcement (Figure [T] left 
inset), and which was soon suggested to be of planetary ori¬ 
gin. Hence, when the event underwent a second much larger 
(factor 10) anomalous brightening 49 days later, OGLE im¬ 
mediately recognized that it would be crucial to fully charac¬ 
terize this second anomaly. In particular, since such violent 
brightening is always caused by the source entering a “caus¬ 
tic” (closed curve of formally infinite magnification - see Fig¬ 
ure 0, a caustic exit was inevitable. Because these caustic 
crossings typically last only an hour or so and can occur at 
any time during the day, OGLE contacted the Microlensing 
Follow Up Network (pFUN) which operates a global network 
of narrow-angle telescopes capable of near 24-hour monitor¬ 
ing of individually important events. With the aid of real-time 
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Fig. 1. — Light curve of the microlensing event OGLE-2013-BLG-0723, including models with and without the planet. Left inset shows the planetary anomaly, 
which includes not just the obvious spike, but also a more extended low level depression. Right inset shows dense coverage of caustic exit, which permitted 
measurement of the Einstein radius Be- 


modeling that accurately predicted the time of the exit, the 
Wise Observatory 0.46m telescope at Mitzpe Ramon, Israel, 
obtained dense coverage of this exit (Figure[Tl right inset). 

The principal data were taken in I band. In addition, OGLE 
obtained some V-band data in order to characterize the mi- 
crolensed source. Finally, /iFUN obtained //-band observa¬ 
tions from the 1.3m SMARTS telescope at Cerro Tololo In- 
teramerican Observatory in Chile. These are not used in the 
present analysis but will be important for the interpretation of 
future follow-up observations to cross-check the lens system 
parameters 

The light curve of OGLE-2013-BLG-0723 shows a system¬ 
atic decline in the baseline. A similar long term linear trend 
(of opposite sign) was seen in OGLE-2013-BLG-0341 and 
was eventually traced to a nearby bright star that was gradu¬ 
ally moving toward (in that case) the source star, so that more 
of its flux was being “captured” in the tapered aperture used 
to estimate the source flux. We searched for such a moving 
bright star by examining the difference of two images, from 
2004 and 2012. We indeed And a dipole from a bright star 
roughly 1.5" away, which is the characteristic signature of 
such moving stars. Having identified the cause of this trend, 
we fit for it and remove it. 

Photometry of the event was done using a customized 
pipeline based on the Difference Imag ing Analysis method 
ilAlard & LuptonlflOOSt lWozniakll200Cii) . Because the photo¬ 
metric errors estimated by an automatic pipeline, cq, are often 


underestimated, the error bars are readjusted by 

e = k(el + elj^^^, ( 1 ) 

where emin is a term used to make the cumulative distribution 
function of as a function of lensing magnification becomes 
linear and the other term k is a scaling factor used to make 
per degree of freedom (dof) becomes unity. The former 
process is needed in order to ensure that the dispersion of data 
points is consistent with error bars of the source brightness 
and the latter process is required to ensure that each data is 
fairly weighted according to error bars. 

3. LIGHT-CURVE ANALYSIS 

The principal features of the light curve are a large double¬ 
horned profile near the peak, which is a signature of binary 
microlenses, and a short-lived spike two months earlier, which 
is characteristic of planetary companions. Modeling these 
features yields the mass ratios and instantaneous locations of 
the three bodies relative to the source trajectory. See Figure|2] 
However, by themselves, these features do not yield either the 
physical or angular scale of the diagram, nor its orientation on 
the sky. Rather, the diagram is scaled to the “Einstein radius”, 
which can be expressed either as an angle 0 e, or as a physical 
length te projected on the observer plane. 

The source star and the Earth’s orbit act as “rulers” to de¬ 
termine these two quantities, as well as the angular orienta¬ 
tion. The source angular radius 0* is known from its observed 
















































UDALSKI ET AL. 


3 


(AU) 


-0.5 0 0.5 1 1.5 



(mas) 


TABLE 1 

Lensing parameters 


Quantity 

Value 

xVdof 

4126.8/4073 

fo (HJD’) 

6484.526 ± 0.037 

Uq 

-0.079 ± 0.002 

fE (days) 

68.48 ± 0.01 

Si 

5.07 ± 0.02 

9i 

3.11 ±0.02 

a (rad) 

-1.195 ±0.003 

•^2 

0.97 ± 0.02 

qi (10-5) 

6.61 ± 0.01 

V’o (rad) 

-4.936 ± 0.005 

p(10-3) 

1.40 ±0.02 

'^E,N 

-0.05 ± 0.01 

fTE.E 

1.35 ±0.02 

ds^/dt (yr“‘) 

0.81 ± 0.02 

dip/dl (yr‘) 

-0.50 ± 0.02 

h 

19.825 ± 0.006 

h 

20.181 ±0.008 


Note. — 1$ and /b represent the baseline magnitudes of the source and 
blend, respectively. 


Fig. 2.— Geometry of the lensing event. The curve with an arrow rep¬ 
resents the source trajectory. Locations of lens components are marked by 
filled dots, where M 3 is the planet, Mj is the planet’s host, and M 2 is the bi¬ 
nary companion. The cuspy closed curves are caustics formed by the lens. 
The inset shows an enlarged view around the caustic that produced the short¬ 
term anomaly. The small empty circle on the source trajectory in the inset 
represents the source size relative to the caustic. We note that the position 
of the planet and the resulting caustic vary in time due to its orbital motion 
around its host. Lower/left axes are expressed by angle on the sky in mas, 
while upper/right axes show projected size of the system in AU. The direc¬ 
tions North and East on the sky are indicated 

temperature and flux. Hence, when it passes over a caustic 
(see Figure |2]), the finite extent of the source smears out what 
would be an extremely sharp peak for a point source (Figure[Tl 
right inset). Second, Earth’s accelerated motion induces dis¬ 
tortions in the light curve, which fix both te relative the size 
of Earth’s orbit (1 AU) and the orientation of the diagram rel¬ 
ative to that of Earth’s orbit, namely the ecliptic. Hence, we 
are able to display both the angular scale and physical orien¬ 
tation in Eigure 2. The mass M and distance Di of the lens 
are then given by 


-L = ^+l. 

Dl rE Ds’ 


4GM 




( 2 ) 


The first equation follows from simple geometry. Once the 
distance is known, the second equation follow from Einstein 
bending angle a = AGM where b is the impact parame¬ 
ter of the lens-sourc e approach. S ee the graphical presentation 
of the formalism in iGouldl ( |2000|) . 


3.1. Microlens Modeling 

Eull modeling of the light curve ultimately requires a to¬ 
tal of 19 parameters, including 14 parameters to describe the 
lensing system and five flux parameters. Some of these man¬ 
ifest themselves in quite subtle effects. The modeling there¬ 
fore proceeded in several distinct phases. Eirst, we determine 
the principal parameters describing the binary-lens system by 
removing the data for 25 days around the planetary pertur¬ 
bation and fitting the remaining light curve to seven system 
parameters plus four flux parameters (two for each observa¬ 
tory). These nine include three parameters to describe the un¬ 
derlying point-lens event (fo,Mo,fE). three to describe the per¬ 
turbations induced by the binary companion (si, g'l, a), one to 


describe the ratio of the source radius to the Einstein radius, 
p = 0*/0E. and the two-dimensional microlens parallax vector 

TTe = (TTE.A'iT^'E.fi)- 

Here fo is the time of closest approach to a reference point 
in the geometry, which is chosen as the “center of magnifica¬ 
tion” of the caustic traversed by the source, mq is the impact 
parameter in units of 0 e. and Ie is the time required for the 
source to cross Be- We choose 0e as the Einstein radius of 
the mass associated with the caustic, which is the total mass 
for the close-binary case, and the mass of the binary compo¬ 
nent located closer to the source trajectory for the wide-binary 
case. This convent ion leads to similar Ie for the two cases, e.g. 
(iGould et al.ll20l4 . 

The three binary parameters are the separation of the com¬ 
ponents in units of Oe, the ratio of the companion to primary 
component masses, and the angle of the source trajectory rel¬ 
ative to the binary axis at fo (since parallax causes this to 
change with time). 

The amplitude of the parallax vector is tte = AU/fe, while 
its direction is that of the lens-source relative proper mo¬ 
tion (in the geocentric frame). In contrast to the lens-source 
trigonometric relative parallax TTrei = tte^e, which is a scalar 
quantity, the microlens parallax requires two parameters be¬ 
cause the lens is not directly observed, so the orientation of 
the lens-sou rce separation re l ative t o the ecliptic is not known 
a priori. See lGould & Horn^(l2Q13h . Normally, one would not 
include parallax in the initial fit, but in this case it is so large 
that the projected Einstein radius is substantially smaller than 
Earth’s orbit, te ~ 0.35 AU, so that there is no even approxi¬ 
mately good solution without including tte in the fit. 

With this parameterization, we find an excellent fit for the 
wide-binary case (s > 1) but only a crude approximation to 
the observed light curve for i < 1. We then incorporate two- 
parameter orbital motion, essentially the instantaneous rela¬ 
tive transverse velocity of the components in Einstein-radius 
units. This refinement results in negligible improvement for 
the i > 1 solution, but dramatic improvement in the i < 1 
solution, so that the two models provide qualitatively simi¬ 
lar approximations of the data. However, it is found that the 
s < 1 case is unlikely due to various reasons as discussed in 
Appendix. 

We then restore the data from the planetary perturbation and 
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Fig. 3.— Disti'ibution of the microlensing parameters of the best-fit solution, whose central values and errors are shown in Table SI. Color coding indicates 
points on the Markov Chain within 1 (red), 2 (yellow), 3 (green), 4 (cyan), 5 (blue) cr of the best fit. 


add five additional parameters to the fit; (i 2 ,^ 2 j' 0 Oj 7 ), where 
7 = {ds 2 /dt/s 2 ,d'ip/dt). Here {s 2 ,q 2 ) are the separation and 
mass ratio of the planet relative to its host (the smaller com¬ 
ponent of the binary), i/jq is the angle between binary axis and 
the planet-host axis at fo, and 7 is the two dimensional pro¬ 
jected velocity of the planet relative to its host. Table [T] shows 
the best ht parameters. In Figure [2 we also present the dis¬ 
tribution of in the parameter space in order to show the 
uncertainties and correlations between the parameters. 

3.2. Origin of Parallax and Orbital Signals 

The microlens parallax signal is quite strong (Ax^ > 1100), 
but since the amplitude tte = 2.7 is one of the largest ever 
detected and hence quite unusual. Therefore, it is prudent to 
check whether this signal is coming mainly from one part of 
the light curve, which would be suspicious, or from the entire 
magnihed portion of the light curve, as one would expect. 

Figure |4] shows that in fact the signal is coming from the 
whole magnihed light curve. It also shows that planetary or¬ 
bital signal, while much weaker (Ax^ = 41), is also broadly 
based. As expected, the orbital signal is more compact than 
the parallax signal. 

3.3. Degeneracy 

Because the light curve is well covered, the parameters are 
well determined. However, the solution is still subject to three 
well known discrete degeneracies. First, there is a degeneracy 
between pairs of companion mass ratio and normalized sepa¬ 
ration {q,s), with one lying inside the Einstein radius (s < 1) 
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Fig. 4.— Cumulative distributions of from the best-fit (orbit+parallax) 
model for the standard (green curve) and parallax (magenta curve) models. 

and the other outside. This degeneracy is fully understood 
at h rst (iGriest & SahzedahllT9^ lDominiklll99^ and second 
(lAnll2005h order for the case that the binary is approximated 
as static. We hnd that within this context there is no “close” 
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Fig. 5.— Location of the source star in the color-magnitude of stars in the 
nearby field. Also marked ai‘e the positions of the centroid of the red clump 
giant stars and of the blend. For the first two, the error bars are smaller than 
the points. For the last, the magnitude is secure but the en'or in the color (not 
shown) is both large and uncertain, being at least 0.2 mag. 

solution, i.e., Ax ^ 2000 between close and wide solutions. 
However, if both close and wide binary models are permitted 
two degrees of orbital motion, i.e. ds^jdt and dil’/dt, then 
the wide solution shows essentially no improvement while the 
close solution improves dramatically and matches the wide 
solution. In Appendix, we discuss the mathematical and phys¬ 
ical nature of this new variant of the close/wide degeneracy 
and show that it implies that the wide solution is strongly fa¬ 
vored. We also present additional corroborating arguments 
that this is the case. We therefore consider only the wide so¬ 
lution. 

Second, if the source lay on the ecliptic, the entire solution 
could be “flipped” relative to the ecliptic without changing 
the light curve. Since the source lies just 4° from the ecliptic, 
this degeneracy is close to exact, with the preferred solution 
favored by only Ax^ = 38. While this is enough to clearly 
favor one solution, the main point is that the two solutions 
yield almost identical physical implications. 

Finally, Figure |2] shows the small planetary caustics just to 
the right of the source trajectory. In principle, there could 
be another solution with these caustics just to the left. How¬ 
ever, this geometry is unable to reproduce the long slow “de¬ 
pression” that follows the planetary anomaly and is due to a 
demagnification “valley” that follows the axis defined by the 
two small caustics. Hence, the solution is both unique and 
measured with good precision. 

3.4. Characteristics of the Source Star and Blended Light 

We determine the angular source radius 0 ^ based on its de - 
reddened color (V — I)o and brightness Iq (IYoo et al.l l2004h . 
which we determine from their offsets from the c entroid of the 
red giant clump r(V-/) ,/1n rinmn = (1.06,14.45) (iBensbv et al.l 
1201 UlNataf et al.ll2013h . under the assumption that the source 
and the clump giants suffer similar extinction. Figure |5] gives 
a graphic representation of this procedure, and also shows the 
position of the star (or ensemble of stars) that are blended 


with the source. That is, in each band, the light-curve flux is 
modeled as F{t) = fsA{t) + ft, where A is the magnification. 
The source and blend magnitudes are then derived from the 
fit parameters /j and ft,, respectively. 

We find (V—I,I)o = (0.52,17.09). From this color and flux, 
we derive = 0.98 ±0.() 8 /xas by converting (V-I) ^ (V- 
K) (iBessell & Bretd 1988h and then app lying a color/surface 
brightness relation (iKervella et al.ll2004l) . The derived color 
is well blueward of the overwhelming majority of Galactic 
bulge stars, except for horizontal branch stars, which are sub¬ 
stantially brighter. Hence, the source is almost certainly in the 
Galactic disk. Disk sources are rare in general (a few percent) 
partly because the source must be farther than the lens with 
the lenses being mostly in the bulge or fairly far in the disk, 
and partly because the column density of bulge stars is higher 
than for disk stars. However, given that this particular lens 
lies far in the foreground {Di ^ 0.5 kpc), it is not unexpected 
for the source to lie in the disk. 

Based on the source color and magnitude, we estimate its 
distance = (7 ± 1) kpc. The exact distance plays no role 
in the determination of 0* because this depends only on the 
flux F oc L/D'^, not the luminosity L and distance separately, 
provided that the source lies behind essentially all the dust. 
Given that the source sits 360 pc above the Galactic plane, 
this is almost certainly the case. The uncertainty in the source 
distance also enters into the uncertainty in the lens distance 
through Equation (O. In practice, however, this leads to an 
error in source distance of only 4%, which is small compared 
to other sources of errors. 

Finally, we note that the source color indicates a temper¬ 
ature T ^ 6100K, which leads to an estimate of the limb- 
darkening co efficient (nee ded for modeling the caustic exit) 
of M = 0.495 (lClaretil200(]h . 

In most microlensing events, the blended light ft, is just a 
nuisance parameter, but in the present case it is potentially 
very important. We first note that while the /-band magnitude 
is secure, the error in the {V - /) color (not shown) is quite 
large. First, the formal error of the fit is 0.2 mag. In addition, 
in contrast to the source flux, which is determined solely from 
changes in brightness, the blend flux depends on the accuracy 
of the baseline photometry. This can be corrupted by bright 
nearby blends, which are plainly visible on the images, and by 
stars that are too faint to be identified. Hence, we consider the 
color to be only indicative and do not make use of it in further 
discussion. However, the I band measurement is too bright to 
be the OGLE-2013-BLG-0723LAby about 1.5 mag. 

To further understand the origin of this light, we measure 
the position of the source relative to the “baseline object”, 
which is an unresolved blend of the source with the star or 
stars giving rise to fb- We find A0base = 55 mas. Since the 
blend contains 43% of the baseline light, this implies that it 
is naively offset from the source by A^biend = A0base/O.43 = 
125 mas. However, while the source position can be mea¬ 
sured extremely well because it becomes highly magnified, 
the “baseline object” is subject to contamination from unseen 
stars as well as the wings of nearby bright stars. We con¬ 
clude that the centroid of blend light is separated from the 
source by A0biend ^ 125 mas. Because the surface density of 
stars I < 20.1 in this field is < 0.1 arcsec”^ (iHoltzman et al.l 
Il998h . the probability that this blended light lies so close to 
the source (and lens) by chance is < 1%. It is therefore most 
likely a companion to either the source or lens. This issue can 
be resolved by obtaining AO images at two epochs to deter- 
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TABLE 2 

Physical parameters 


Quantity 

Value 

Mass of the planet 

Mass of the host 

Mass of the binary companion 

Distance to the lens 

Projected planet-host separation 

Projected separation between binary components 
Relative lens-source proper motion (geocentric) 
Relative lens-source proper motion (heliocentric) 

0.69 ± 0.06 M® 
0.031 ±0.003 Mq 
0.097 ± 0.009 Mq 
0.49 ± 0.04 kpc 

0.34 ± 0.03 AU 
1.74±0.I5 AU 

3.75 ±0.33 mas/yr 
14.5 ±1.3 mas/yr 


mine whether its proper motion is consistent with the source 
or lens. If it proves to be associated with the lens, then this 
is a four-body system consisting of a brown dwarf orbited by 
a Venus-mass planet, with both orbiting a very low-mass star, 
and with this entire group orbiting another star that is roughly 
double the mass of the star that contributed to the microlens- 
ing event. 

3.5. Physical Parameters 

The principal physical parameters are given in Table |2l 
From the determined masses of the individual lens compo¬ 
nents, it is found that the lens is a triple system that is com¬ 
posed of a Venus-mass planet (0.69 M^) orbiting a brown- 
dwarf host (0.031 Mq) that is accompanied by a low-mass 
star (0.097 Mq). The projected planet-host separation is 
0.34 AU, while the separation between the binary components 
is 1.7 AU. The system is located at a distance 490 pc toward 
the Galactic center. 

Note that we give the proper motion in both the geocentric 
frame in which the measurements were actually made, and the 
heliocentric frame, which will be useful for comparison with 
future observations. Since is known, the conversion from 
geocentric to heliocentric frames is straightforward. 

3.6. Which Binary Component is Host to the Planet? 

From Figure|2]or Table[Tl the planet’s projected position is 
5.07/0.97 = 5.2 times closer than the stellar companion. If the 
planet were orbiting the more massive companion, this would 
occur with probability < 4%. In addition, in order to ensure 
stability of the 3-body system, the brown dwarf would have to 
be behind (or in front of) the star-planet system by a at least 
factor of 3 larger distance than the projected separation of the 
brown dwarf and star. This chance alignment further reduces 
the probability by a factor 3^. Hence, while this configura¬ 
tion is possible in principle, it has extremely low probability. 
A very similar argument applies to a circumbinary configura¬ 
tion, in which the planet’s orbit is at least three times larger 
than the separation between the brown dwarf and star, but just 
happens to be seen in projection close to the brown dwarf. 

4. DISCUSSION 

OGLE-2013-BLG-0723LBb is a missing link between 
planets and moons. That is, its host OGLE-2013-BLG- 
0723LB, being a brown dwarf in orbit around a self-luminous 
star, is intermediate between stars and planets, in both size and 
hierarchical position. Moreover, the scaled mass and host- 
companion separation of OGLE-2013-BLG-0723LB/Bb are 
very similar to both planets and moons in the solar system. 
See Table [3l Both Uranus and Callisto are believed to have 
formed in the cold outer regions of their respective accretion 
disks, and are mostly composed of the raw materials of such 


regions; ice with some rock. In the case of Uranus, it is be¬ 
lieved to have been formed closer to the current location of 
Saturn (10 AU) and to have migrated outward. In the table, 
the companion-host separations are scaled to the host mass. 
This is appropriate because the “snow line”, the inner radius 
at which icy solids can form (2.7 AU in the solar system) in¬ 
creases with host mass, probably roughly linearly. A plausi¬ 
ble inference from the first three lines of Table 2 is that these 
processes scale all the way from solar-type stars hosting plan¬ 
ets, to brown dwarfs hosting “moon/planets”, to giant planets 
hosting moons. 

However, while OGLE-2013-BLG-0723LBb is almost cer¬ 
tainly orbiting OGLE-2013-BLG-0723LB today, we cannot 
guarantee that it was born there. Planets in relatively close 
binaries (like OGLE-2013-BLG-0723LA,B) can jump from 
one star (or in this case, brown dwarf) to the other, if the 
system is dynamically perturbed (iKratter & PeretsI [201211 . If 
this is what happened to OGLE-2013-BLG-0723LBb, then 
would have been captured from an orbit around OGLE-2013- 
BLG-0723LA to its current orbit around OGLE-2013-BLG- 
0723LB. It would then be similar to Triton, which was cap¬ 
tured by Neptune out of solar orbit. Table [3] shows that the 
Neptune-Triton system in its current configuration is also not 
terribly different than a scaled version of OGLE-2013-BLG- 
0723LB/Bb. 

It will be interesting to model the dynamics of the triple sys¬ 
tem OGLE-2013-BLG-0723LA,B,Bb, but it would be prema¬ 
ture to do so. This is because there is likely a fourth member 
of this system that is more massive and luminous than the 
other three, and which is separated by of order 100 AU or 
less. If confirmed, this member could significantly impact the 
system dynamics via the Lidov-Kozai mechanism. The prin¬ 
cipal evidence for this fourth member is that there is excess 
light in the aperture that is too bright to be the primary star 
OGLE-2013-BLG-0723LA, and is so closely aligned with the 
source that it is unlikely to be a random interloper. Because 
the source flux is known from the //-band light curve (which 
is not included in the present analysis), it will be straightfor¬ 
ward to subtract out the source light from a high-resolution 
adaptive optics (AO) image and determine the exact location 
of the blended light. If well separated, then its association 
with the lens can be established from its proper motion, de¬ 
rived from observations at two epochs, which can be com¬ 
pared with Table 121 If directly superposed in the AO images, 
then it is almost certainly associated with either the source or 
the lens because of the extremely low probability of a random 
interloper. The two possibilities can be distinguished based 
on the apparent color because the source is highly reddened 
by dust, but the lens sits in front of the great majority of the 
dust. 

OGLE-2013-BLG-0723LBb is the second terrestrial mass 
planet orbiting one member of a binary system composed of 
two very low mass objec ts. The other was O GLE-2013-BLG- 
0341LBb. Gould et al. (iGould et al.ll2014|) showed that if all 
low-mass stars were in such systems, roughly one should have 
been detected. They noted that no strong statistical inference 
could be derived from this fact because, for example, if only 
1 % of low-mass stars were in such systems, there would be 
a 1% chance to find one. Because the question about the fre¬ 
quency of these systems was being framed a posteriori, a 1% 
effect is not particularly significant. 

However, OGLE-2013-BLG-0723LBb reproduces three of 
the four elements that went into their estimate of the fre¬ 
quency of such detections, namely small probability of pass- 
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TABLE 3 

Companion/Host pairs 


Host/Companion 

^comp /^host 

(lO'A 

^c-h / ^host 

(AU/M©) 

Uranus/Sun 

4.4 

19 

OGLE-2013-BLG-0723LBb/B 

6.6 

11 

Callisto/Jupiter 

5.7 

13 

Inton/Neptune 

21 

46 


ing by the planetary caustic, small probability of crossing the 
central (binary) caustic, and location in a high-cadence OGLE 
fields. The fourth criterion, source magnitude brighter than 
Is = 18.5, is not satisfied, since Is ^ 19.8. If the limit were ex¬ 
tended to Is < 20, the total number of available events would 
roughly double. Hence, the prediction would rise to two. With 


two predicted and two detected, the probability that these are 
extremely rare drops roughly as the square, i.e., roughly 10“"^ 
probability of this occurring by chance if only 1% of these 
stars were in such systems. Hence, while there is still no com¬ 
pelling evidence that these systems are extremely common, at 
least they are not extremely rare. Eurther, the current detec¬ 
tions remain consistent with the hypothesis that these systems 
are ubiquitous. 
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acknowledge support from NSE AST-1103471. B.S.G., 

A.G., and R.W.P. acknowledge support from NASA grant 
NNX12AB99G. 


APPENDIX 

RESOLUTION OE THE CLOSE/WIDE DEGENERACY 


In the main text of the paper, we stated that the wide solution was consistent with a static binary while the close solution 
required substantial orbital motion to match the wide-solution light curve. At first sight this seems qualitatively reasonable since 
close binaries have shorter periods than wide ones. However, only one of these solutions coiTesponds to the physical system, 
while the other is merely a mathematical reflection of it. Hence, the real question is not whether each solution is physically 
plausible but what is the relative likelihood where each solution, if it were real, could generate the other as a doppelganger. 

Before addressing this, some comments must be in order as why this degeneracy is appearing in this case and how likely it 
will be to appear generically. In the simplest cases two solutions appear because to first order, the quadrupole expansion of the 
close-bi nary potential and the tidal expansion of t he wide-binary potential are identical and so (to this order) produce identical 
caustics jGriest & Safizedahll998t[Dominikll999h . Contemporaneously, it was noticed that in more typical cases, the degeneracy 
could remain ex tremely severe even though the caustics were not actually identical and, indeed, appeared rotated relative t o each 
other (Eig. 8 of (I Afonso et al.ll200^ L This was ultimately explained by carrying out the analysis to second order dAnl200-5l) . That 
is, an additional degree of freedom (basically, rotation) could compensate for the slower convergence of the potential. Eor wide 
binaries, the potential converges more slowly as the mass ratio of the companion increases and as its separation gets closer. This 
is well illustrated by small mass-ratio, e.g., planet/star, systems. Topologically, the two caustics are the same, but geometrically 
they are very different. The planet has a relatively weak effect on the stellar potential, meaning that the star’s caustic (“central 
caustic”) is well-described by low-order terms and hence subject to the close/wide degeneracy. By contrasts, the planetary caustic 
is strongly perturbed and not subject to this degeneracy (unless the separation is very large). 

In the present case, the caustic feature in the light curve of the wide solution is associated with the lower-mass companion, and 
the mass ratio is relatively large. Hence, one expects that (as in the planetary case) higher ord er terms will play a relatively large 
role. And indeed they prove to be too large to permit the degeneracy analyzed by (lAnll2005l) . Orbital motion then provides the 
additional degree of freedom to compensate for the non-trivial higher-order terms. 

Consider then a wide binary that is qualitatively similar to the present case. Regardless of details, its light curve can always 
be mimicked by including the appropriate amount of orbital motion. On the other hand, consider a (real) close binary with the 
(s,q) parameters predicted by the wide binary but with so far unspecified orbital parameters. Can this system “generate” a static 
wide binary with the same light curve? The answer is: only if its orbital parameters happen to lie in a naiTow range around those 
predicted by the wide pair. Otherwise, the static wide binary would have been consistent with a wide range of close-binary orbital 
motions. Because Ax^ is quite large, the fine-tuning requirements of the close solution are severe. This is the principal reason 
that we discount it. 

However, in addition, if the close solution is spurious and merely tuned to the mathematical requirements of matching the wide 
solution, it is likely to have characteristics that are atypical of real systems, which can be a further indic ation of its spurious 
nature. We note two such features. Eirst, the ratio of projected kinetic to potential energy (iDong et al.ll2009h is quite low 


UpotA M/Mg 


(Al) 


compared to typical values 0.1 ^ P ^ 0.5. The three physical effects that can produce low /3 are 1) it is near apocenter of an 
extremely eccentric orbit, 2) our line of sight happens to be almost tangent to the orbital velocity vector, or 3) we happen to view 
the system edge on, so that its semi-major axis is much larger than its projected separation. Since the unprojected ratio typically 
has values Ekin/Epot 0.5 reducing this number by a factor ^ 20 via projection effects requires some fine tuning. However, in 
the present case there is an additional complication: the third option generates the need for an additional low probability event. 
That is, if the two stars are extended along the line of sight, then either this projection must be very great to allow the planet to 
have a stable orbit around one star (actually brown dwarf), to the planet itself must be seen in projection next to the pair despite 


















the fact that its orbit is at least three times that of the extended binary. Finally, we note that the source star is four times fainter 
in the close model. Since its color (so spectral type) is still the same, this places it twice as far as in the wide solution, where the 
density of such stars is about a factor 40 smaller. This is somewhat compensated by the eight times larger phase space, but still 
does further reduce the probability of this solution by a factor 5. 

All the arguments given above make the close-binary solution very unlikely but they do not exclude it. However, this solution 
can also be tested by direct imaging. Since there are compelling reasons to conduct such imaging, possibly at several epochs, no 
special effort is required to do this. If these observations reveal a blend that is substantially brighter than the source, then the close 
solution will be confirmed. The expected baseline magnitudes of the source and blend are (/s,/b) = (21.28,19.42), respectively. 
If the wide solution is correct, we expect the blend and source to be comparable (depending on their /-// colors). 
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